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Abstract — This paper presents a new framework for multime- 
dia streaming that integrates the application and network layer 
functionalities to meet such stringent application requirements 
as delay and loss. The coordination between these two layers 
provides more robust media transmission even under severe 
network conditions. In this framework, a multiple description 
source coder is used to produce multiple independently-decodable 
streams that are routed over partially link-disjoint (non-shared) 
paths to combat bursty packet losses. We model multi-path 
streaming and propose a multi-path selection method that chooses 
a set of paths maximizing the overall quality at the client. Overlay 
infrastructure is then used to achieve multi-path routing over 
these selected paths. The simulation results show that the average 
peak signal-to-noise ratio (PSNR) improves by up to 8.1 dB, if the 
same source video is routed over intelligently selected multiple 
paths instead of the shortest path or maximally link-disjoint 
paths. In addition to PSNR improvement in quality, the end- 
user experiences a more continual streaming quality. 

I. Introduction 

With the advances in audio-visual encoding standards 
and broadband access networks, multimedia communications 
(MMC) has become closer to reality. The popularity of the 
Internet further stimulated the demand for MMC services. 
However, the underlying infrastructure, i.e., routing and trans- 
port protocols, and control and signalling mechanisms, has 
not been sufficiently developed to support such challenging 
applications. Therefore, early efforts largely failed, and this 
experience encouraged researchers to seek new techniques 
to support delay-sensitive multimedia applications over lossy 
networks. 

In this paper, we propose a novel framework for mul- 
timedia streaming over packet networks. While generating 
several streams with a multiple description (MD) source coder, 
we transmit these streams over diverse paths by using an 
overlay infrastructure. Current routing algorithms are mainly 
destination-based and they use the hop counts as their distance 
metric. Although, this idea performs well for data applications, 
it generally underperforms for real-time multimedia applica- 
tions. One of the key mechanisms needed for incorporating 
MDC into routing algorithms is the ability to efficiently select 
multiple paths between the server and client. An easy solution 
would be to use totally link-disjoint paths. However, these 
paths or the shortest paths do not necessarily make the best use 



of error-resiliency properties of MDC. Moreover, totally link- 
disjoint paths are rarely available between two end hosts. The 
routing algorithms for MDC-enabled multimedia streaming 
system should choose paths in such a way that the overall 
streaming quality is maximized or equivalently, the overall 
distortion is minimized. 

The contemporaneous work about using MDC with path 
diversity is [1]. Based on this work, Apostolopoulos et al. 
discuss MD streaming within content delivery networks in [2]. 
The authors study video streaming using different reference 
picture selection schemes and path diversity in [3], [4] and 
show that the quality can be increased in doing so. Also, the 
study [5] investigates the performance of five different traffic 
dispersion techniques for high bandwidth MPEG-2 streaming 
data. All of these research studies propose methods to increase 
the reliability of the streamed video. However, to emulate 
multi-path transport they assume the availability of multiple 
statically independent paths between the end hosts, which is 
not feasible in many cases. Furthermore, they neither provide 
a particular network model nor a path selection method. In 
this study, we address these open problems and propose a 
method to select optimal multiple paths as a function of 
the network parameters, media characteristics and application 
requirements. 1 

The rest of the paper is organized as follows. In Section II, 
we give an overview of MDC. Section HI introduces the 
envisioned network model in this work. Multi-path selection 
is discussed in Section IV. Simulation results are presented in 
Section V. Finally, Section VI concludes the paper. 

II. Multiple Description Coding Theory 

Multiple Description Coding (MDC) is a source coding 
that generates several descriptions such that different levels of 
reconstruction qualities can be obtained from different subsets 
of these descriptions. In contrast to Layered Coding, there is 
no hierarchy among the descriptions so that each description 
may be independently decoded. This property makes MDC 
highly suitable for packet networks where no prioritization 

1 It should be mentioned that the proposed multi-path selection method is 
only optimal within the context of the network models and objective functions 
that are considered in this paper. 
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Fig. 1. A typical multiple description codec with two descriptions. 

exists among the packets. MDC provides this robustness at the 
expense of a small reduction in the compression efficiency. 

A typical MDC system is shown in Fig. 1. A sequence of 
source symbols {Xk}^ =l is input to the multiple description 
(MD) encoder. The MD encoder produces two descriptions, 
which are then transmitted over possibly different channels. 
If both descriptions are received intact at the client, then 
we use the central decoder to produce the highest quality 
signal. However, if only one description is received error-free, 
then, the corresponding side decoder is used. We represent the 
reconstructed signal at the i th decoder by {X^}^^ and the 
distortion by Ri denotes the rate in terms of the number 
of bits per source sample on the i th channel. 

Our objective is to estimate the end-to-end video quality in 
terms of path parameters. We start our derivation by relating 
the rate of the descriptions to the reconstructed signal quality, 
or equivalently, the distortion. For the sake of simplicity, we 
consider only two descriptions generated at the source. We 
define the average distortion at the i th decoder as 




i = 0,1 and 2, (1) 



where d(.) can be any distortion measure. In this paper, we 
take d{.) as the squared Euclidian norm. 

To develop analytic expressions for the average distortion, 
we need to choose a media type and source model. In the 
rest of the paper, we assume that the source is a compressed 
video stream. However, note that methodology presented here 
applies to any source with its specific distortion model. That is, 
the proposed path selection method chooses the paths resulting 
in the maximal quality irrespective of the type of the streamed 
media and its distortion model. 

In compressed-video applications, the distribution of Dis- 
crete Cosine Transform (DCT) coefficients can be modelled 
as Laplacian. For a Laplacian source with variance cr 2 , the 
rate-distortion function is upper and lower bounded by [6] 

V2- 2 * ^ Du^euM) < ° 2 2- 2Ri • (2) 

7T 

In Section IV, we will use (2) to derive a relation between the 
source rate (Ri) and streaming distortion. 

There are various techniques for generating multiple de- 
scriptions. One of the most straightforward methods is time- 
domain partitioning. For instance, separating the even and 
odd numbered frames in a video sequence into two groups 
produces two descriptions. Likewise, the sequence may also 
be partitioned in the spatial-domain. Other popular approaches 



are Multiple Description Scalar Quantization (MDSQ) and 
Multiple Description Transform Coding (MDTC). In MDSQ, 
multiple quantizers are used to generate descriptions whereas 
correlating transforms are used in MDTC. For brevity, we will 
not go into further details of MDC. For an excellent survey, 
interested readers are referred to [7]. 

III. Envisioned End-to-end Network Model 

In this study, we assume the availability of an overlay 
network that spans several nodes including the server and 
client similar to the Resilient Overlay Networks (RONs) and 
X-Bone [8], [9]. In such an overlay network, we are able to 
achieve MD streaming by application- layer routing. Although 
it is not as efficient as network-layer routing, for overlays, 
application-layer routing does not require any network support, 
hence, can easily be implemented in today's Internet. Here, 
the overlay nodes, called O-nodes, are logically connected to 
each other to build an overlay network that will be exploited 
to enable multi-path routing. With this aspect, we first find the 
optimal set of paths from the server to the client. Then the path 
information, i.e., addresses of the O-nodes to be traversed, is 
encoded into the packet pay load. This information is used by 
the O-nodes to forward the incoming packet to the next O- 
node on the path. Between the O-nodes, the packets are routed 
by the underlying routing mechanism such as the shortest path 
routing. Hence, the overall procedure is transparent to the 
physical network. 

Besides enabling multi-path routing, each O-node contin- 
uously monitors the statistics of the flows between other O- 
nodes and itself. To discover the overlay topology and link 
characteristics, every O-node exchanges information among 
other O-nodes by a routing protocol. The throughput, packet 
loss rate, delay and jitter statistics are then used to compute 
application-specific routes. The details of implementing such 
protocols are further discussed in [8]. 

A. Definitions 

Now, we will define the network and the link parameters. 
Recall that, we are operating on an overlay network. Hence, 
a link stands for a logical link. 

Definition 1: A network is defined by the tuple (N,L) 
where N = {Ni, N 2i N n } represents the nodes in the 
network and L is an n x n matrix. L(u, v) = 1 if a direct 
link, denoted by l UiV , exists between nodes N u and N v and 
L(u,v) = 0, otherwise. 

Definition 2: A path Vs,c is defined by a set of nodes and 
links connecting nodes N$ and Nc- The set of the nodes on 
this path is represented by Ns,c Q N and the set of the links 
on this path is denoted by £s,c- 

Let *Psc anc * c b e anv tw0 P atns - These paths are said 
to be totally link-disjoint if and only if C l s c and £| c do not 
have any common element. If the sets C\ c and C 2 S c intersect, 
the links belonging to both paths are referred as joint links, 
whereas the rest of the links are referred as the disjoint links. 



B. Link Parameters 

The parameters for the link l UtV are defined as follows: 

• &«,v denotes the bandwidth between nodes N u and N v . 

• Pu, v denotes the packet loss probability between nodes 
N u and N v . p u , v values for different links are assumed 
to be independent. 

• tu lV denotes the minimum delay between nodes N u and 
N v . This includes processing, transmission and propaga- 
tion delays, but not varying queuing delays. 

• ju,v denotes the average jitter between nodes N u and N v . 
This jitter is mainly due to the varying queuing delays. 
Together with t UtVi j UjV determines the forward trip time. 

IV. Multi-path Selection 

In this section, we introduce the multi-path selection method 
step by step. Let's consider two paths, Vg C and V\ c , 
between the nodes Ns and Nc, which are referred as a path 
pair. We do not impose that the paths should be different, 

Vsc mav oe me same as ^sc ® ur S oa ^ 1S t0 estimate the 
end-user quality in terms of the link parameters making up 
these paths. In other words, we are looking for a function / 
that determines the average distortion at the client, D c , based 
on the link parameters. That is, 

D c = f{K y v,Pu,vitu,viju,v) for V/ UjV e £s,c> * = 1)2. (3) 

Recall that our ultimate goal is to find the path pair that 
minimizes the cost function, /. For two descriptions, we 
have four possible cases determined by the description on- 
time arrivals. Let Z) 0| o (i\i) denote the distortion when 
both descriptions arrive intact on time (are lost or delayed). 
Similarly, let £>o,i (^i,o) denote the distortion when the first 
(second) description arrives intact on time, but the other one 
is lost or delayed. Given these distortions, we can write the 
average distortion at the client as 

r Al a 

D c = Prob {Both received on time} x JD 0 ,o (4) 

/ A s 

+ Prob {1 st received ontime&; 2 nd lost or delayed} xD 0 ,i 
+ Prob {1 st lost or delayed & 2 nd received on time} x D x 0 

N v ' 

+ Prob {Both lost or delayed} x£> M , 

s V ' 

where the success probabilities, Fo,o — Pi,u will be derived 
later in this section. The shared parameters due to the joint 
links play a critical role in the path pair selection process. 
Hence, the selection process should consider all possible path 
pairs. One cannot evaluate the paths individually, rank them, 
and then choose the best two. This does not necessarily give 
the optimal pair due to the dependencies between the paths. 
For the sake of simplicity, throughout the section we drop the 
subscripts and superscript and use i to denote the i th path, 
where i € {1, 2}. We start our derivations with distortion. 
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Fig. 2. State transition diagram for the Gilbert-Elliott model. 

A. Bandwidth - Distortion Relation 

Let Bi denote the end-to-end available bandwidth on the 
path Vi. Clearly, B< is the minimum bandwidth over all the 
links on V%. That is, B{ = min{b UjV : V/ UjV € £*}• In 
Section II, we have given the upper and lower bounds for the 
distortion of a Laplacian source. Now, we will give the relation 
between Ri in (2) and B{. Assuming that the i th description is 
encoded with a rate of Ri bits per pixel (bpp), and the source 
video has a resolution of W x H (pixels/frame) at a frame 
rate of F (frames/second), we have 

K' WxfxFxc (bpp) ' (5) 
where c is a known constant that depends on the chroma 
subsampling format. Typically, a streamed video is in 4:2:0 
format and the corresponding c is 1.5. Since there is only a 
coefficient difference between the lower and upper bounds in 
(2), for our purposes, it suffices to give the distortion as a 
function of the aggregated bandwidth of the path pair (B) as 

D(B) = <t 2 2~ 2x wxAfxc . (6) 

Given this relation, we can compute the distortion terms in 
(4) by plugging the corresponding B values into (6). First, 
for Do,o the aggregated bandwidth is typically the sum of 
the individual bandwidths, i.e., B = B\ + B<z provided that 
the joint links have available bandwidths higher than this 
sum. However, if this is not the case, B is then limited by 
the bottleneck joint-link capacity. Then, to compute £) 0 ,i and 
£>i,o» we take B = f?i and B = B 2 , respectively. Finally, 
Di,! is computed by letting B = 0. 

B. Computation of Success Probabilities 

The next step in the computation of (4) is to find the success 
probabilities (Po,o — Pi t i)- Each of these probabilities are 
comprised of two components. The first component computes 
the probability that the description arrives at the receiver in a 
finite duration, i.e., it is not lost during the transmission. We 
will refer this component as arrival probability. On the other 
hand, for the second component we are interested in comput- 
ing the probability that the description arrival occurs before 
a pre-specified deadline. This component will be referred as 
on-time delivery probability. 

B.l. Arrival Probability 

We adopted a well-known two-state Markovian Gilbert- 
Elliott (GE) model [10], [1 1] to describe the temporal behavior 
of packet losses on a link. The state transition diagram of the 
GE model is given in Fig. 2. In this figure, G (good) and B 
(bad) represent the packet arrival and loss states, respectively 
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Fig. 3. Path segregation method not only reduces the topology complexity 
but also allows us to compute the arrival probabilities accurately. 

and the corresponding steady-state probabilities are given by 
ttq — and 7r B = -^^0. For the link / u>v , we compute the 
transition probabilities by 

&u,v = ~rz and Pu,v = 7 — , (7) 

■t Pu t v ■*- / u i v 

where is the average burst length on this link. Typically, 
L UtV is a number between 4 and 7 for applications running over 
UDR It is important to note that if temporal dependency of 
bursty packet losses was ignored, the optimal solution would 
be to send all descriptions over the best single path. 

The GE model provides a simple expression for the charac- 
terization of individual links. However, derivation of the end- 
to-end arrival probabilities on partially link-disjoint paths is 
not trivial. Next, we focus on this crucial step. 

B.l.a. Path Segregation 

The GE model successfully captures bursty loss character- 
istics of a flow. However, when several flows are interleaved 
on the same link, the burst observed by a particular flow 
may become less significant. As enough number of foreign 
packets are introduced in between two consecutive packets of 
a flow, packet loss events become less and less correlated and 
eventually, the GE model reduces to a Bernoulli model. 

Consider the topology depicted in Fig. 3 where two can- 
didate paths split at node N z , and merge later at node Nq. 
The joint links between the server and the segregation point 
(N z ) form a sub-network denoted by Network^*. Beyond 
the node N z , the disjoint links traversed by each path until 
the node Nq are also grouped to form the corresponding sub- 
networks that are denoted by Network^ isjoint . 

Assume that the source generates two correlated descrip- 
tions, X\ l and X l 2 \ with a period of At where U denotes 
the time at which the descriptions are generated. The value of 
At is chosen large enough so that the descriptions generated 
at the time instants U and U 4- At will not fall into the 
same burst period, i.e., the loss events for X[ l and x{ i+At 
will be independent of each other. Further assume that each 
description is one packet and transmitted as soon as it is ready. 
Consequently, the packets belonging to X{* and X% will arrive 
at node Nz approximately at the same time (back-to-back). 
However, once the paths split at node N z , these two packets 
will be routed over different paths. The chances are that these 
packets will not arrive at the node Nq within a typical burst 
duration. Moreover, beyond the node Nq, it is not likely that 
these packets will again meet in the same burst duration. 
From this observation, we deduce that once the paths split, 



the loss events on the rest of the paths should be considered 
independent despite a possible merge at a later node. This 
suggests us to employ a joint GE model for Networkj 0in£ , 
and individual Bernoulli models for Network^ oin£ . We refer 
this method as path segregation. 

We have simulated several topologies to check the validity 
of path segregation. We compared the packet loss burst size 
histograms generated by the original and approximated topolo- 
gies. As expected as the delay difference between the paths 
became comparable to the average burst duration, the model 
error reduced almost to zero. In a real network, the disjoint 
paths will always have enough delay difference since they will 
experience uneven delay and jitter. 

B.l.b. Link Aggregation 

Generally, the joint sub-path in Network^* is a combina- 
tion of several Markovian links. Although it is straightforward 
to compute the arrival probabilities by considering each of 
these links individually, this computation will be greatly sim- 
plified if we can approximate the combination of these links 
by a single Markovian link. However, for the approximated 
link one has to find the joint GE parameters, namely aj 0 i nt 
and (3 joint- Recall that consecutive links' packet loss rates are 
independent. Hence, we can directly compute the packet loss 
rate for Network Join£ as 

P joint = 1- IJ (l-p Uy v)- (8) 

Vi u ,t>€ Network j 0 i nt 

However, estimating the joint average burst length, Lj 0int 
is not as trivial. This is because bursty periods in a stream 
accumulate as the stream goes over more links, i.e., effective 
average burst length increases. Although there are several 
approaches to handle this problem [10], [11], they usually 
require higher order state analysis that renders the proposed 
approach impractical. For this reason, we estimate Lj oint with 
the following heuristic: 

ljoint — x ^ ] Pu,v x L UjV . (9) 

P joint yf lutv€ Network ioiflt 

Once we have the Pj 0 int and Lj oint , we compute aj 0int and 
Pjoint by plugging p joint and L joint into (7). This method 
is referred as link aggregation since it combines several 
Markovian links into a single one. The simulation results show 
that the heuristic in (9) successfully estimates Lj 0int . 

On the other hand, since the disjoint sub-paths in 
Netwovk l dis j oint are modelled by Bernoulli, we can directly 
convert these links into a single Bernoulli link whose loss 
probability is given by 

Pdisjoint = l- II U-Pu.v). (10) 

VJu.veNetwork^.^ 

Finally, the arrival probability of a path is the product of the 
arrival probabilities on its joint and disjoint sub-paths. 



B.2. On-time Delivery Probability 

The second step in computing the success probabilities is 
to derive on-time delivery probabilities. Let T, denote the 
minimum end-to-end delay for the path V*. We can write Ti as 
Ti = Sviu v ed Obviously, we compel T { to be smaller 
than the pre-specified deadline period. Otherwise, it means that 
all descriptions transmitted over this path will be obsolete by 
the time they are received. 

Let Ji denote the average end-to-end jitter on path V{. 
We divide J, into two components, jj oint and j disjpinV as 
the jitters experienced in Network^ 0 i n t and Network^ oinl , 
respectively. Assuming that jitters on the consecutive links are 
independent, we compute Jj as follows: 

Ji= 3u,v+ 3 U > V ' 

Vi UiU €Networkj 0 i„t Vf u Network* . 

s ' > v < 

Jjoint jji 

Jdiajoint 

Let TdeadUne represent the maximum tolerable delay for a 
target application. Generally, different paths have different 
end-to-end delays, and hence, different tolerances for jitter. 
We define i^ ax as the maximum jitter that can be tolerated 
on path Vi. By definition, tj nax = T dead u ne - Then, it is 
easy to see that the on-time delivery probability can also be 
given as the probability that Ji is smaller than t l max . 

Finally, we have all the information required to write 
success probabilities (P 0 ,o - Pi t i) in terms of path parameters. 
Recall that, in (4), P 0 ,o computes the probability that both 
descriptions arrive at the client intact on time. To find this 
value, we first compute the arrival probability of transmitting 
both descriptions intact over Network JO i n( and individual 
arrival probabilities of transmitting each description intact over 
Network^ SJoint . Since these events are independent, we mul- 
tiply them to find the arrival probability of both descriptions. 
Secondly, we compute the probability that the descriptions will 
arrive at the client before the pre-specified deadline. Finally, 
we multiply this value with the arrival probability to get Po,o- 

On the other hand, P 0 ,i considers the cases where the first 
description is received intact on time and the second descrip- 
tion is (i) lost in Networkj 0 i nt or in Network^ oinf or (ii) 
delayed beyond the play-out time although transmitted intact. 
Hence, we separately compute the corresponding probabilities 
of these three events and then add them up to get P 0 ,i. Pi,o 
can also be computed in the same manner. Lastly, P% y i is the 
complement of other given probabilities. The corresponding 
equations are given in (12) - (15) where we denote the joint 
GE parameters by a and (5 instead of aj 0 i nt and ft joint, 
respectively. 

V. Simulation Results 

We have conducted several simulations to demonstrate the 
efficacy of the proposed multi-path selection method for a 
real-time streaming application. To this end, we generated 
a random Internet topology and by varying the topology 
parameters we compared different streaming methods in terms 
of video quality. 



Po.o =(Orc(l " a) 2 + n B 0{l - a)) X (1 - p diajoint ) x (1 - p 2 diajoint )) 02) 

X PrOb{jjoint + jdiajoint ^ iox Midpoint + 3 disjoint < *max) 

Po.i =((t G (1 - 0)0 + *b0*) x (1 - p l diajoint )) (13) 

X Prob(j joint + j dia j oint < 'max) 

+ ((tt g (1 - a) 2 + * B 0{1 - a)) x (1 - P disjoint ) X p 2 diajoint ) 

X PrOb(j joint +3 dia joint < *max) 

+ ((tt g (1 - a) 2 + n B p(l - a)) x (1 - p di3joint ) x (1 - p 2 disjoint )) 

X Pr0b(j join t + disjoint < 'max ^ Jjoint + j dia j 0 int > <max) 

Pi,o =((7r G a/3 + ttb(1 - 0)0) x (1 - P 2 diajoint )) (14) 

X Prob{j joint + ^disjoint < 'max) 

+ (0tg(1 - a) 2 + v B 0(l - o)) x p l diBjoint x (1 - p disjoint )) 

X Prob(J joint + J'L joint < t ma x) 

+ ((tt g (1 - a) 2 + 7T B 0(l - a)) x (1 - p diajoint ) x (1 - p 3 diajoint j) 

X Pr0b{j joint +3 dia joint > 'max Md jjoint + Jdisjoint < 'max) 

Pi,i =1 - Po,o - Po,i - Pi.o (15) 

We used two standard test sequences Table Tennis and 
Flower Garden in our simulations to stream from a server 
to a client with a delay tolerance of 200 ms. Each sequence 
comprised of 150 frames. The resolutions were 352 x 240 pix- 
els and the frame rate was 30 frames per second. To make 
quantitative comparisons, we used the PSNR measure given 
by PSNR = 10 x log 10 (^§^), where MSE stands for the 
mean squared error. Among various multiple description video 
coders, we preferred the time-domain partitioning method with 
two descriptions. 

In the simulations, the single description (SD) and multiple 
description (MD) encoded streams were produced with a 
standard MPEG-2 encoder at the available rates. Note that we 
were not interested in any particular video encoding scheme 
or concealment technique that could be applied to increase the 
streaming quality. Here, we performed a simple concealment 
technique that repeated the information from the last available 
frame in case of a bursty loss. For the isolated single errors, we 
concealed the missing frame with its neighboring frames (the 
other description in the MD case). Naturally, more intelligent 
encoding and concealment techniques would result in higher 
streaming qualities for both SD and MD encoded streams. 

A. Internet Topology Simulations 

In order to generate a random Internet topology, we used 
GT-ITM [12]. The resulting topology consisted of 192 nodes 
with 4 mesh-connected transit domains, 8 nodes per transit 
domain, 1 stub per transit node and 5 nodes in a stub 
domain. 2 An overlay node (O-node) was associated with each 
transit domain node to enable multi-path routing. Notice that 

2 Since we were interested in streaming from a stub domain to another one 
(where the server and client were located, respectively), we kept the number 
of stub domains per transit domain node small. 




Fig. 4. An Internet topology with 4 transit domains. 



multi-path routing capability was intended within only transit 
domains, not stub domains. An illustrative topology with 4 
transit domains is shown in Fig. 4. 

For the link capacities, we assigned 1.0 Gbps to transit- 
to-transit edges; 100 Mbps to transit-to-stub edges; 5 Mbps 
and 2 Mbps to stub-to-stub edges. The available bandwidth 
on each link (b UjV ) was chosen randomly in the range of 
20%-80% of its capacity. The packet loss probabilities (p UyV ) 
were randomly assigned proportional to the link utilizations 
in the range of 0.1% - 1.0% for all links except one link in 
each domain (which was also randomly chosen to emulate 
congestion), whose loss rate was between 10%-20%. On the 
other hand, we assigned delays (t u ,v) in tne range of 10 - 
30 ms for the transit-domain links and 5 - 10 ms for the 
stub-domain links. Lastly, average jitter values for all links 
(ju,v) were randomly chosen a value between 0-10 ms. These 
averages were then used to introduce exponentially distributed 
jitter for each packet. 

On this topology, we randomly chose a server-client pair 
such that the server and client were separated by at least 
two transit domains. We compared three different streaming 
methods namely, the shortest path, maximally link-disjoint 
path and optimal multi-path streaming. For the shortest path, 
we chose the one with minimum number of hops and if 
this metric was equal for two or more paths, we selected 
the one with the minimum delay. On the other hand, for 
the maximally link-disjoint path pair, we first generated and 
evaluated all possible path pairs. Then, we looked for the pairs 
with minimum number of overlapping links to minimize the 
statistical dependency between the paths and then chose the 
pair with minimum total end-to-end delay. Notice that in doing 
so, we not only approximated the totally link-disjoint pairs, 
which were supposed to exist in [3], [4], but also selected 
shorter ones to avoid quality degradation due to the delayed 
packets. Lastly, for optimal multi-path streaming, based on the 
network statistics we selected the pair among all possible ones 
that minimized (4). While SD encoded video was used in the 
first method, we used the MD encoded video in the second 
and third methods. 

We ran 10 different simulations for each test sequence by 



choosing different congested links at each run. 3 This analysis 
allowed us to observe the variations in the performances of 
three streaming methods. In the last simulation (Run #10) 
for the TABLE Tennis sequence, we deliberately congested 
the link by which the stub domain including the server was 
connected to its transit domain. In this case, as expected, all 
three streaming methods performed poor since all streams had 
no choice other than to be routed over the same congested 
link. Obviously, path diversity becomes less effective when 
the congestion occurs on the joint links before the multiple 
paths split. The results for the two test sequences are tabulated 
in Tables I and II, respectively. In each row, we present 
the average PSNR of the streamed video and total number 
of lost and delayed (obsolete) packets for a different run. 
Evidently, optimal multi-path streaming outperforms shortest 
path and maximally link-disjoint path streaming in all of these 
independent runs. There are prominent conclusions that can be 
drawn from these results. We elaborate on them next. 

TABLE I 

Simulation results for the Table Tennis sequence 





Shortest Path 


Max. Link-Dis. Path 


Optimal Multi-Path 


Run# 


PSNR 


Lost/Del. 


PSNR 


Lost/Del. 


PSNR 


Lost/Del. 


1 


33.1 


4/ 1 


33.2 


4/2 


33.5 


3/2 


2 


25.7 


28/3 


30.2 


15/3 


33.0 


6/5 


3 


32.5 


6/ 1 


33.2 


5/2 


33.2 


3/4 


4 


32.4 


7 / 1 


33.3 


5/ 1 


33.5 


3/2 


5 


25.7 


28/3 


30.1 


18/3 


33.2 


3/4 


6 


29.3 


18/2 


30.1 


19/4 


33.0 


5/5 


7 


31.6 


8/2 


33.5 


4/ 1 


33.5 


4/ 1 


8 


25.2 


32/3 


30.2 


17/2 


33.3 


3/3 


9 


32.4 


6/2 


32.0 


10/3 


33.3 


3/3 


10 


25.8 


28/2 


27.2 


25/4 


27.3 


24/5 


Avg. 


29.4 


16.5 / 2.0 


31.3 


12.2/2.5 


32.7 


5.7 / 3.4 



TABLE H 

Simulation results for the Flower Garden sequence 





Shortest Path 


Max. Link-Dis. Path 


Optimal Multi-Path 


Run# 


PSNR 


Lost/Del. 


PSNR 


Lost/Del. 


PSNR 


Lost/Del. 


1 


18.8 


30/3 


21.9 


16/2 


24.1 


4/3 


2 


19.1 


27/3 


21.7 


17/3 


24.1 


3/4 


3 


18.9 


29/2 


21.6 


18/2 


23.9 


4/4 


4 


24.0 


5/ 1 


24.1 


5/2 


24.3 


3/3 


5 


23.6 


6/2 


24.4 


4/2 


24.4 


3/3 


6 


18.9 


28/3 


21.8 


15/3 


23.8 


6/5 


7 


23.8 


6/2 


23.9 


8/3 


24.4 


3/3 


8 


23.9 


6/ 1 


24.3 


5/ 1 


24.5 


4/2 


9 


20.3 


20/2 


21.0 


20/4 


23.8 


5/4 


10 


23.6 


6/2 


23.7 


5/3 


23.8 


4/4 


Avg. 


21.5 


16.3 / 2.1 


22.9 


11.3/2.5 


24.1 


3.9 / 3.5 



B. Discussion 

When we carefully examine the results, we observe that 
the shortest path streaming often degrades the quality to an 
unacceptable level. This occurs when a link on the shortest 
path is heavily congested. Naturally, the maximally link- 
disjoint path streaming is able to reduce the adverse effect of 
the congested link by exploiting the path diversity and improve 
the average quality over the shortest path streaming. However, 

3 We repeated the simulations 100 times to obtain statistically reliable 
results. Also, we chose an average burst length, L, of 4. 



Table Tennis - Simulation #5 
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Fig. 5. The PSNR comparisons of the shortest path and optimal multi-path 
streaming methods for Table Tennis. 

Flower Garden - Simulation #2 
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Fig. 6. The PSNR comparisons of the shortest path and optimal multi-j 
streaming methods for Flower Garden. 



its performance is still below the optimal multi-path streaming. 
Hence, we conclude that minimizing the dependency between 
the paths without considering their characteristics does not 
necessarily result in the maximal streaming quality. 

Another interesting point is the distribution of individual 
frame qualities. The plots in Fig. 5 and Fig. 6 show the PSNR 
values of each frame for the videos streamed over the shortest 
path and optimal path pair. These plots are extracted from the 
5 th and 2 nd simulations for the Table TENNIS and FLOWER 
Garden sequences, respectively. As seen from the plots, the 
poor video quality persists for the SD encoded stream during 
a bursty loss. Moreover, the wide fluctuations further degrades 
the perceptual quality. However, the video quality slightly 
deviates from its average value for the MD encoded stream 
unless both descriptions are lost simultaneously. This feature 
provides smoother streaming experience to the client. 

The results also show that optimal multi-path streaming 
experienced a slightly higher number of delayed packets 
compared to other methods. This is mainly because of the 
fact that in this method, a longer but less congested link may 
be preferred over a shorter but more congested link as long as 
the delay tolerance allows to do so. On the other hand, when 
the shortest path (the maximally link-disjoint paths) is not 



congested, Le. 9 it is a good path, the proposed method may also 
choose this shortest path (the maximally link-disjoint paths) as 
well. Because of this flexibility, optimal multi-path streaming 
performs better than its opponents. Lastly, we observed the 
tendency of multi-path routing in load balancing. Rather than 
overloading the shortest path, the total load is shared between 
different paths. 

VI. Conclusions 

In this paper, we modelled multi-path streaming and pro- 
posed an optimal multi-path selection method for MD stream- 
ing. This method yields the maximal video quality at the 
client by evaluating the network conditions and application 
requirements. The simulations show that the optimal multi- 
path selection attains considerably large quality improvements 
over both the shortest and maximally link-disjoint path stream- 
ing. At the implementation side, we achieve the MD streaming 
by the use of overlay infrastructure. In doing so, we do not 
advocate any particular congestion or routing mechanism in 
the underlying network. However, additional improvement can 
be achieved if the network nodes become multi-path routing 
aware, e.g., one of the descriptions may be prioritized over the 
others to prevent concurrent description losses. 

In this paper, we paid attention to the two-path case due to 
its tractability and practicability, however it is straightforward 
to generalize the proposed model to more than two paths at 
the expense of increased computational complexity. Also, the 
framework can be adapted to other streaming media such as 
audio and 3-D graphics. 
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